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(57) Abstract 

In a system and method of digital wireless communi- 
cation between a base station (B) and a subscriber unit (S), 
a spatial channel characterized by a channel matrix (H) cou- 
ples an adaptive array of (MT) antenna elements at the base 
station (B) with an adaptive array of antenna elements (MR) 
at the subscriber station (S). The method comprises the use 
of spatio-temporal coding (TRELLIS ENCODER), training 
symbols (TRAINING SYMBOL INJECTION), and frequency 
domain deinterleaving (INTERLEAVER), At the receiver, a 
matched de-interleaver (DE-INTERLEAVER) transforms the 
space-frequency sequence back into a serial signal stream. A 
maximum likelihood detector (VECTOR MAXIMUM LIKE- 
LIHOOD DETECTOR) generates the recovered information 
stream. 
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High Capacity Wireless Communication 
Using Spatio-Temporal Coding 

RELATED APPLICATIONS 

This application claims priority from U.S. provisional 
applications 60/025,227 and 60/025,228, both filed 08/29/96. 
Both applications are hereby incorporated by reference. 

FIELD OF THE INVENTION 

This invention relates generally to digital wireless 
communication systems. More particularly, it relates to using 
antenna arrays by both a base station and a subscriber to 
significantly increase the capacity of wireless communication 
systems . 

BACKGROUND OF THE INVENTION 

Due to the increasing demand for wireless communication, it 
has become necessary to develop techniques for more 
efficiently using the allocated frequency bands, i.e. 
increasing the capacity to communicate information within a 
limited available bandwidth. This increased capacity can be 
used to enhance system performance by increasing the number of 
information channels, by increasing the channel information 
rates and/or by increasing the channel reliability. 

FIG. 1 shows a conventional low capacity wireless 
communication system. Information is transmitted from a base 
station B to subscribers Si, . . . ,S 9 by broadcasting 
omnidirectional signals on one of several predetermined 
frequency channels. Similarly, the subscribers transmit 
information back to the base station by broadcasting similar 
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signals on one of the frequency channels. In this system, 
multiple users independently access the system through the 
division of the frequency band into distinct subband frequency 
channels. This technique is known as frequency division 
5 multiple access (FDMA) . 

A standard technique used by commercial wireless phone systems 
to increasing capacity is to divide the service region into 
spatial cells, as shown in FIG . 2. Instead of using just one 

10 base station to serve all users in the region, a collection of 
base stations Bi , ...,67 are used to independently service 
separate spatial cells. In such a cellular system, multiple 
users can reuse the same frequency channel without interfering 
with each other, provided they access the system from 

15 different spatial cells. The cellular concept, therefore, is 
a simple type of spatial division multiple access (SDMA) . 

In the case of digital communication, additional techniques 
can be used to increase capacity. A few well known examples 

20 are time division multiple access (TDMA) and code division 
multiple access (CDMA) . TDMA allows several users to share a 
single frequency channel by assigning their data to distinct 
time slots. CDMA is normally a spread-spectrum technique that 
does not limit individual signals to narrow frequency channels 

25 but spreads them throughout the frequency spectrum of the 
entire band. Signals sharing the band are distinguished by 
assigning them different orthogonal digital code sequences. 
These techniques use digital coding to make more efficient use 
of the available spectrum. 

30 

Wireless systems may also use combinations of the above 
techniques to increase capacity, e.g. FDMA /CDMA and TDMA/ CDMA . 
Although these and other known techniques increase the 
capacity of wireless communication systems, there is still a 
35 need to further increase system performance. Recently, 
considerable attention has focused on ways to increasing 
capacity by further exploiting the spatial domain. 
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One well-known SDMA technique is to provide the base station 
with a set of independently controlled directional antennas, 
thereby dividing the cell into separate sectors, each 
controlled by a separate antenna. As a result, the frequency 
reuse in the system can be increased and/ or cochannel 
interference can be reduced. Instead of independently 
controlled directional antennas, this technique can also be 
implemented with a coherently controlled antenna array, as 
shown in FIG. 3. Using a signal processor to control the 
relative phases of the signals applied to the antenna 
elements, predetermined beams can be formed in the directions 
of the separate sectors. Similar signal processing can be 
used to selectively receive signals only from within the 
distinct sectors. 



In an environment containing a significant number of 
reflectors (such as buildings), a signal will often follow 
multiple paths. Because multipath reflections alter the 
signal directions, the cell space experiences angular mixing 
and can not be sharply divided into distinct sectors. 
Multipath can therefore cause cochannel interference between 
sectors, reducing the benefit of sectoring the cell. In 
addition, because the separate parts of such a multipath 
signal can arrive with different phases that destructively 
interfere, multipath can result in unpredictable signal 
fading. 



In order to avoid the above problems with multipath, more 
sophisticated SDMA techniques have been proposed. For 
example, U.S. Pat. No. 5,471,647 and U.S. Pat. No. 5,634,199, 
both to Gerlach et al . , and U.S. Pat. No. 5,592,490 to Barratt 
et al. disclose wireless communication systems that increase 
performance by exploiting the spatial domain. In the 
downlink, the base station determines the spatial channel of 
each subscriber and uses this channel information to 
adaptively control its antenna array to form customized beams. 
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as shown in FIG. 4A . These beams transmit an information 
signal x over multiple paths so that the signal x arrives to 
the subscriber with maximum strength. The beams can also be 
selected to direct nulls to other subscribers so that 
5 cochannel interference is reduced. In the uplink, as shown in 
FIG. 4B, the base station uses the channel information to 
spatially filter the received signals so that the transmitted 
signal x' is received with maximum sensitivity and 
distinguished from the signals transmitted by other 
10 subscribers. In this approach the same information signal 
follows several paths, providing increased spatial redundancy. 

In the uplink, there are well known signal processing 
techniques for estimating the spatial channel from the signals 

15 received at the base station antenna array, e.g. by using a 
priori spatial or temporal structures present in the signal, 
or by blind adaptive estimation. If the uplink and downlink 
frequencies are the same, then the spatial channel for the 
downlink is directly related to the spatial channel for the 

20 uplink, and the base can use the known uplink channel 
information to perform transmit beamforming in the downlink. 
Because the spatial channel is frequency dependent and the 
uplink and downlink frequencies are often different, the base 
does not always have sufficient information to derive the 

25 downlink spatial channel information. One technique for 
obtaining downlink channel information is for the subscriber 
to periodically transmit test signals to the base on the 
downlink frequency rather than the uplink frequency. Another 
technique is for the base to transmit test signals and for the 

30 subscriber to feedback channel information to the base. If 
the spatial channel is quickly changing due to the relative 
movement of the base, the subscriber and/ or reflectors in the 
environment, then the spatial channel must be updated 
frequently, placing a heavy demand on the system. One method 

35 to reduce the required feedback rates is to track only the 
subspace spanned by the time-averaged channel vector, rather 
than the instantaneous channel vector. Even with this 
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reduction, however, the required feedback rates are still a 
large fraction of the signal information rate. 

Although these adaptive beamforming techniques require 
5 substantial signal processing and/or large feedback rates to 
determine the spat ial channel in real time, these techniques 
have the advantage that they can navigate the complex spatial 
environment and avoid, to some extent, the problems introduced 
by multipath reflections. As a result, an increase in 

10 performance is enjoyed by adaptive antenna array systems, due 
to their use of the spatial dimension. Note, however, that 
while the base station antenna array can make efficient use of 
the spatial dimension by selectively directing the downlink 
signal to the subscriber S, the uplink signal in these systems 

15 is spatially inefficient. Typically, the subscriber is 
equipped with only a single antenna that radiates signal 
energy in all directions, potentially causing cochannel 
interference. These communication systems, therefore, do not 
make optimal use of the spatial dimension to increase 

20 capacity. 

OBJECTS AND ADVANTAGES OF THE INVENTION 

Accordingly, it is a primary object of the present invention 
to provide a communication system that significantly increases 

25 the capacity and performance of wireless communication systems 
by taking maximum advantage of the spatial domain. Another 
object of the invention is to provide computationally 
efficient coding techniques that make optimal use of the 
spatial dimensions of the channel. In particular, it is an 

30 object of the present invention to provide coding techniques 
specially adapted for the case of rapidly fading channels 
where channel state information (CSI) at the transmitter is 
unknown. These and other objects and advantages will become 
apparent from the following description and associated 

3 5 drawings. 
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SUMMARY OF THE INVENTION 

These objects and advantages are attained by a method of 
digital wireless communication that takes maximal advantage of 
spatial channel dimensions between a base station and a 
5 subscriber unit to increase system capacity and performance. 
Surprisingly, the techniques of the present invention provide 
an increased information capacity in multipath environments. 
In contrast, known techniques suffer in the presence of 
multipath and do not exploit multipath to directly increase 

10 system capacity. In brief, the present invention teaches a 
method of wireless communication using antenna arrays at both 
the base and subscriber units to transmit distinct information 
signals over different spatial channels in parallel, thereby 
multiplying the capacity between the base and the subscriber. 

15 In particular, the present invention teaches specific spatio- 
temporal coding techniques that make optimal use of these 
additional spatial subchannels in the case of unknown 
transmitter channel state information. 



20 Generally, the present invention provides a method of digital 
wireless communication between a base station and a subscriber 
unit in the case where channel state information is not known 
by the transmitter. For this purpose a spatio-temporal coding 
structure that exploits the spatial subchannel capacity is 

25 used. In particular, a matrix orthogonal frequency division 
multiplexing (MOFDM) scheme and a space- frequency trellis 
coding system is used at the transmitter, and a space- 
frequency maximum likelihood detector with a channel estimator 
are used at the receiver. With this relatively simple 

30 structure, a MIMO system according to the present invention is 
able to provide a channel capacity several times greater than 
can be achieved in a conventional wireless system using OFDM. 
The inventors also propose an efficient channel estimation 
algorithm for the time varying MIMO channel. 
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DESCRIPTION OF THE FIGURES 

FIG. 1 shows a low capacity wireless communication system well 

known in the prior art. 
FIG. 2 illustrates a known technique of spatially dividing a 

service region into cells in order to increase system 

capacity. 

FIG. 3 illustrates the use of beamforming with an antenna 
array to divide a cell into angular sectors, as is known 
in the art . 

FIGS. 4A and 4B illustrate state-of-the-art techniques using 
adaptive antenna arrays for downlink and uplink 
beamforming, respectively . 

FIGS. 5A and 5B show the parallel transmission of distinct 
information signals using spatial subchannels in downlink 
and uplink, respectively, as taught by the present 
invention . 

FIGS. 6A and 6B are physical and schematic representations, 
respectively, of a communication channel for a system 
with multiple transmitting antennas and multiple 
receiving antennas, according to the present invention. 

FIGS. 7A and 7B are block diagrams of the system architecture 
for communicating information over a multiple- input- 
multiple-output spatial channel according to the present 
invention. 

DETAILED DESCRIPTION 

Although the following detailed description contains many 
specifics for the purposes of illustration, anyone of ordinary 
skill in the art will appreciate that many variations and 
alterations to the following details are within the scope of 
the invention. Accordingly, the following preferred 
embodiment of the invention is set forth without any loss of 
generality to, and without imposing limitations upon, the 
claimed invention. 

As discussed above in relation to FIGS. 4A and 4B, prior art 
wireless systems employing an adaptive antenna array at the 
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base station are multiple-input-single-output (MISO) systems, 
i.e. the channel from the base to the subscriber is 
characterized by multiple inputs at the transmitting antenna 
array and a single output at the receiving subscriber antenna. 
5 Because these MISO systems can exploit some of the spatial 
channel, they have an increased capacity as compared to 
single-input-single-output (SISO) systems that are discussed 
above in relation to FIGS. 1 and 2. It should be noted that 
although the MISO systems disclosed in the prior art provide 

10 an increase in overall system capacity by spatially isolating 
separate subscribers from each other, these systems do not 
provide an increase in the capacity of information transmitted 
from the base to a single subscriber, or vice versa. As shown 
in FIGS. 4A and 4B, only one information signal is transmitted 

15 between the base and subscriber in both downlink and uplink of 
a MISO system. Even in the case where the subscriber is 
provided with an antenna array, the prior art suggests only 
that this capability would further reduce cochannel 
interference. Although the overall system capacity could be 

20 increased, this would not increase the capacity between the 
base and a single subscriber. 

The present invention, in contrast, is a multiple-input- 
multiple-output (MIMO) wireless communication system that is 

25 distinguished by the fact that it increases the capacity of 
both uplink and downlink transmissions between a base and a 
subscriber through a novel use of additional spatial channel 
dimensions. The present inventors have recognized the 
possibility of exploiting multiple parallel spatial 

30 subchannels between a base station and a subscriber, thereby 
making use of additional spatial dimensions to increase the 
capacity of wireless communication. Surprisingly, this 
technique provides an increased information capacity and 
performance in multipath environments, a result that is in 

35 striking contrast with conventional wisdom. 
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FIGS. 5A and 5B illustrate a MIMO wireless communication 
system according to the present invention. As shown in FIG. 
5A, a base station B uses adaptive antenna arrays and spatial 
processing to transmit distinct downlink signals xi , yi 2 , X3 
through separate spatial subchannels to a subscriber unit S 
which uses an adaptive array and spatial processing to receive 
the separate signals. In a similar manner, the subscriber S 
uses an adaptive array to transmit distinct uplink signals 
x'i, x' 2 , x* 3 to the base B over the same spatial subchannels. 
As the multipath in the environment increases, the channel 
acquires a richer spatial structure that allows more 
subchannels to be used for increased capacity. 

It is important to note that the simple assignment of the 
distinct signals to the distinct spatial paths in a one-to-one 
correspondence, as illustrated above, is only one possible way 
to exploit the additional capacity provided by the spatial 
subchannel structure. For example, coding techniques can be 
used to mix the signal information among the various paths. 
In addition, the present inventors have developed techniques 
for coupling these additional spatial dimensions to available 
temporal and/or frequency dimensions prior to transmission. 
Although such coupled spatio-temporal coding techniques are 
more subtle than direct spatial coding alone, they provide 
better system performance, as will be described in detail 
below. 

It is also important to note that the transmit beamforming at 
the base requires knowledge of the downlink channel state 
information. (Similarly, the transmit beamforming at the 
subscriber requires knowledge of the uplink channel state 
information. Because the system is symmetric with respect to 
the base and subscriber, it suffices to discuss one case.) 
Although downlink channel state information can be fed back to 
the base from the subscriber, if the channel is rapidly 
changing, then the demand on the channel capacity to provide 
real time channel information and the demand on the signal 
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processing may make it impractical to implement the system 
under the assumption that transmit channel state information 
is available. Accordingly, the inventors have developed an 
MOFDM coding technique to take advantage of the added spatial 
5 subchannels even in the case of unknown transmitter channel 
state information . 

In order to facilitate an understanding of the present 
invention and enable those skilled in the art to practice it, 

10 the following description includes a teaching of the general 
principles of the invention, as well as implementation 
details. First we develop a compact model for understanding 
frequency dispersive, spatially selective wireless MIMO 
channels in the case where the channels are time invariant, 

15 and then generalize to the case where the channels vary with 
time. We then discuss the theoretical information capacity 
limits of these channels, and propose spatio-temporal coding 
structures that exploit the spatial subchannel capacity in the 
case of unknown channel state information. In particular, a 

20 matrix orthogonal frequency division multiplexing (MOFDM) 
scheme is described. In a preferred embodiment a space- 
frequency trellis coding system is located at the transmitter, 
and a space- frequency maximum likelihood detector with a 
channel estimator are located at the receiver. With this 

25 relatively simple structure, a MIMO system according to the 
present invention is able to provide a channel capacity 
several times greater than can be achieved in a conventional 
wireless system using OFDM. The inventors also propose an 
efficient channel estimation algorithm for the time varying 

3 0 MIMO channel. 

In its preferred implementations, the present invention makes 
use of many techniques and devices well known in the art of 
adaptive antenna arrays systems and associated digital 
35 beamforming signal processing. These techniques and devices 
are described in detail in U.S. Pat. No. 5,471,647 and U.S. 
Pat. No. 5,634,199, both to Gerlach et al . , and U.S. Pat. No. 
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5,592,490 to Barratt et al . , which are all incorporated herein 
by reference. In addition, a comprehensive treatment of the 
present state of the art is given by John Livita and Titus 
Kwok-Yeung Lo in Digital Beamforming in Wireless 
Communications (Artech House Publishers, 1996). Accordingly, 
the following detailed description focuses upon the specific 
signal processing techniques which are required to enable 
those skilled in the art to practice the present invention. 

Consider first a time-invariant communication channel for a 
system with M T transmitting antennas at a base B and M R 
receiving antennas at a subscriber S, as illustrated in FIGS. 
6A and 6B. The channel input at a sample time k can be 
represented by an M T dimensional column vector 

Z(k) = [Zi(k) Z MT (k)]T, 

and the channel output and noise for sample k can be 
represented, respectively, by M R dimensional column vectors 

x(k) = [xi (k) , . . . ,x MR (k) ]T, 
and 

n(k) = [ni(k) n MR (k) 

The communication over the channel H may then be expressed as 
a vector equation 

x(k) = Hz(k) + n(k) , 

where the MIMO channel matrix is 
f hi, l ■ • - hi,M T \ 



H = 



^h MR ,i . . h MR/MT> 



Each matrix element hj.j represents the SISO channel between 
the i th receiver antenna and the j th transmitter antenna. Due 
to the multipath structure of the spatial channel, orthogonal 
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spatial subchannels can be determined by calculating the 
independent modes (e.g. eigenvectors) of the channel matrix H. 
These spatial subchannels can then be used to transmit 
independent signals and increase the capacity of the 
5 communication link between the base B and the subscriber S. 

In the case where the channel matrix H is not fixed in time, 
but changes, it should be represented as a time-dependent 
matrix, H(k) . Moreover, because the multipath introduces time 
10 delays into the various propagation paths, a spatial 
decomposition of H independent of time will result in temporal 
mixing of the signals. It is more appropriate, therefore, to 
perform a more general spatio-temporal analysis of the 
channel . 

15 

Let {zj(n)} be a digital symbol sequence to be transmitted 
from the j th antenna element, g(t) a pulse shaping function 
impulse response, and T the symbol period. Then the signal 
applied to the j th antenna element at time t is given by 

20 

sj (t) = £ Zj (n)g(t-nT) 
n 

The pulse shaping function is typically the convolution of two 
separate filters, one at the transmitter and one at the 

25 receiver. The optimum receiver filter is a matched filter. In 
practice, the pulse shape is windowed resulting in a finite 
duration impulse response. We assume synchronous complex 
baseband sampling with symbol period T. We define no and (v+1) 
to be the maximum lag and length over all paths 1 for the 

30 windowed pulse function sequences {g(nT - Tj)}. To simplify 
notation, it is assumed that no = 0, and the discrete- time 
notation g(nT - %i) = gi (n) is adopted. 

When a block of N data symbols are transmitted, N+v non-zero 

35 output samples result. Denoting k as the block index for the 
k th channel usage, k(N+v) is the discrete time index for the 
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first received sample, and (k+l)(N+v)-l is the time index for 
the last received sample. The composite channel output can now 
be written as an M R - (N+v) dimensional column vector with all 
time samples for a given receive antenna appearing in order so 
that 



x(k) = [xi (k{N+v) ) , . . . ,xi ( (k+1) (N+V) -1) # 

x MR (k(N+V) > , . . . ,x Mr ( (k+1) (N+v) -1) ]T # 



with an identical stacking for 
Similarly, the channel input 
vector written as 



the output noise samples n(k) . 
is an M T -N dimensional column 



z(k) = [zi (k(N+V) ),..., zi (k(N+V)+N-l) 

z MT {k(N+V) ) z MT (k(N+v)+N-l) ]T / 



The spatio-temporal communication over the channel H(k) may 
then be expressed as a vector equation 

x(k) = H(k)z(k) + n(k) , 



where the MIMO time -dependent channel matrix 
/Hi # i(k> . . . H i/MT (k) \ 



H(k) = 



H MR ,i(k) * • H MR , MT (k) 



is composed of SXSO sub-blocks Hjj (k) . 

To clearly illustrate the effect of multipath, the channel can 
be written as the sum over multipath components 



H(k) = X 



a R ,i<e R# ,)l 

a R,M R 



^(k)G r [a T/1 (G T/f )I ... a T/MT (6 T|f )l]. 
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In this equation, aR,j(0R,f) is the gain response of the j th 
receiver array element due to angle of arrival Qr,i of the f th 
multipath signal, aT,i(0T ( f) is the gain response of the i th 
transmitter array element due to angle of departure Qt,1 of the 
fth raultipath signal. Kif(k) is the diagonal time varying 
channel fading parameter matrix given by 

«i(k) = diag [ Pr(k(N+v)) < (k+1) (N+v) -1) ], 



10 



and the Toeplitz pulse shaping matrix G{ is given by 



g,(0) 0 0 0 . .. 0 

g,(v) ... g ? (0) 0 0 
0 g ? (v) ... g ? (0) 0 0 



Ql (V) 



g,(0) 



0 0 g ? (v) 



We will now discuss the information capacity for the spatio- 
temporal channel developed above. The following analysis 
assumes that the noise n(k) is additive white Gaussian noise 
(AWGN) with covariance a 2 I v +i - Each channel use consists of an 

N symbol burst transmission and the total average power 
radiated from all antennas and all time samples is constrained 
to less than a constant. 

Write the singular value decomposition (SVD) of the channel 
matrix as H (k) =V H (k) A H (k) u£ (k) , with the n th singular value 
denoted XH r n(k)- Write the spatio-temporal covariance matrix 
25 for z (k) for block index k as R z (k) with eigenvalue 
decomposition R z (k) =V Z (k) A z (k) Uz (k) , and eigenvalues X Z/n (k). 



15 



20 



It can be demonstrated that, if the case where the 

instantaneous channel state information is known at both the 

30 transmitter and receiver, then the information capacity for 

14 
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the time-varying discrete- time spatio-temporal communication 
channel defined above is given by 



C = E 



( 

| l og fx + ^.n(k)|X Hrn (lc)| 2 ^ 



K 

I 

n=l 



V 



where X, z , n (k) is given by the spatio-temporal water-filling 
solution, E(-) is the expectation operator, and K is the 
number of finite amplitude singular values in H(k) . 

For the case where only the receiver has instantaneous channel 
state information, it is not possible to adapt the transmitter 
for each block. Nevertheless, it is possible to find the time 
invariant transmitter covariance which maximizes the capacity 
for the worst case channel possibilities. For any given 
transmitted signal covariance matrix R z , the worst case 
channel would place all of the time average energy in the rank 
1 subspace defined by the smallest eigendirection in R 2 . This 
game theoretic problem leads to a spatially uncorrelated 
transmitter covariance solution R 2 = g x Mt , where p T is the 

maximum average block transmission power. This transmitter 
covariance is used for completely unknown point to point 
channels and broadcast channels. For this case of unknown CSI 
at the transmitter, it can be shown that a white space- time 
transmission distribution gives a channel capacity 



r 



C = E 



K 

I log 

n=l 



V 



1 + 



PTl^H.n(k) P 

M T a 2 



J 



By analyzing the ranks of the matrices in the path 
decomposition of the time varying channel H(t), it can be 
demonstrated that the maximum number of finite amplitude 
parallel spatio-temporal channel dimensions, K, that can be 
created to communicate over the far field time-varying channel 
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defined above is equal to min{ N-L, (N+v) -M R , N-M T }, where L 
is the number of multipath components. Thus, multipath is an 
advantage in far-field MIMO channels. If the multipath is 
large ( L » 1 ) , the capacity can be multiplied by adding 
5 antennas to both sides of the radio link. This capacity 
improvement occurs with no penalty in average radiated power 
or frequency bandwidth because the number of parallel channel 
dimensions is increased. In practice, an adaptive antenna 
array base station, such as that described by Barratt et al., 

10 is modified to implement a coding scheme, as described below, 
which exploits these additional dimensions. In particular, a 
signal processor is designed to perform a spatio-temporal 
transform of information signals in accordance with the above 
equations so that they may be transmitted through the 

15 independent parallel subchannels and decoded by the 
subscriber . 

In constant or slowly varying channels, it is often possible 
to send training sequences to the receiver and communicate 

20 channel state information (CSI) back to the transmitter in a 
manner that accurately tracks time variations. In such cases, 
the transmitter can implement a coding solution which 
approaches the theoretical capacity limits. The MIMO 
communication problem becomes more difficult when the channel 

25 fades rapidly in time as is the case with portable wireless 
communication in the microwave frequency bands . It then 
becomes impractical to feed back CSI f rom the receiver to the 
transmitter due to the information bandwidth required to 
update the channel state in real time. It is highly desirable 

30 in such cases to have a channel coding technique that exploits 
the spatial dimension of the MIMO problem without requiring 
any CSI at the transmitter. Such a coding technique has been 
devised by the present inventors and is described in detail 
below. 

35 

Given the time-varying channel defined by H(t), it is 
theoretically possible to create a coding system consisting of 
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a spatio-temporal encoder, and a spatio-temporal maximum 
likelihood decoder. The obvious difficulty with such a system 
is the complexity of the decoder. The complexity of the 
spatio-temporal decoder can be greatly reduced, however, by 
using a matrix orthogonal frequency division multiplexing 
(MOFDM) structure according to the present invention. The 
complexity reduction occurs because inter-symbol interference 
(ISI) is eliminated from each OFDM sub-channel. 

The MOFDM channel structure is derived under the assumption 
that the channel is block time invariant over a block of N+2v 
symbol periods. Under this assumption, the channel fading 
matrix R f (k) can be replaced by the scalar fading variable 
(3?(k). Note that the block time invariant assumption is 
reasonable provided that the block duration (N+2v)T « Ap, 
where Ap is the correlation interval for the channel fading 
variable. (The correlation interval is defined here as the 
time period required for the fading parameter time- 
autocorrelation function to decrese to some fraction of the 
zero-shift value.) 

For MOFDM, N data symbols are transmitted during each channel 
usage. However, a cyclic prefix is added to the data so that 
the last v data symbols form a preample to the N data symbol 
message block. By discarding the first and last v data 
symbols at the receiver and retaining only N time samples at 
the channel output, the new MIMO channel H(k) has a block 
cyclic structure: 



L 



a R,l< 9 RW> X 



H(k) = £ft(k) 



[a T/1 (9 T j )I . . . 




a R,M R (Or ,? > x 



where the cyclic pulse shaping matrix G(k) is given by 
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O, = 



g,(0> 0 ... 0 

g f (V - 1) ... g,(0) 0 
g, (V) g f (V - 1) ... g,<0) 



g f (v) 



g f (v) g f (v - l) 



... g f (D 

0 g, ( v ) 
0 

• . g,(0) 



10 



The MOFDM channel model can now be derived as follows. First 
post multiply H(k) with the N-M T x N-M T block diagonal inverse 
discrete Fourier transform (IDFT) matrix f* (M T ) w here each 
diagonal block is the unitary N x N IDFT matrix F* . The next 
step is to premultiply by a similar N-M R x N-M R block diagonal 
DFT matrix F< m r> where the diagonal submatrices F are N x N 
DFT matrices. Pre- and post-multiplication by permutation 
matrices Pr and then gives the decomposition of the channel 
into discrete discrete Fourier transform (DFT) frequency 
domain sub-channels M n (k), as follows: 



15 



9t(k) = V (N)P R F< M R)H(k)F* ( M T>P T 

fXi(k) 0 > 
o "w N {k)> 



Each channel 5t n (k} is independent of the other frequency 
domain sub-channels. Just as in the case of scalar OFDM, the 
20 cyclic prefix allows the large time domain channel to be 
decomposed into many smaller parallel frequency domain 
channels. The received vector signal 1X> n (k) for each frequency 
domain spatial sub-channel can then be expressed as 

25 

tt n (k) = W n (k)Z n (k) + Tln(k) . 

where Z n (k) is the subchannel transmitted signal and TL n (k) is 
the subchannel noise. A system architecture implementing this 
30 channel structure is shown in FIG. 7A. 
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The spatial sub-channels can also be expressed as 

L 

K„(k) = X P(<k) 9?,n «r,i aj ( , 



where 9,, n is the DFT of the sequence {g, (k) } evaluated at DFT 
index n. At each frequency index, the DMMT channel is due to a 
weighted sum over L rank-1 outer products of the frequency- 
invariant receive and transmit array response vectors. The 
weighting is determined by the frequency invariant path fading 
values and the Fourier transform of the delayed pulse shaping 
function. This reveals a highly structured nature for the time 
varying space- frequency channel spectrum. 

In the case of rapidly fading channels where CSX is not 
available at the transmitter, the appropriate transmitter 
distribution is a spatially and temporally white transmitter 
sequence. Nevertheless, as seen from the above channel 
decomposition, the use of cyclic signal structures allows the 
determination of a channel structure that can still be 
exploited to improve capacity. Therefore, a practical 
subchannel coding method which approximates a white 
distribution is desired. Although many variations are 
possible, the following description is focused on a 
particularly simple strategy involving a one dimensional 
trellis coding structure. 

In MOFDM, a space- frequency code is transmitted. Given M T 
transmitting antenna elements and the MOFDM subchannel 
decomposition of % (k) , a codeword sequence c^) of constraint 
length N C D can be viewed as q spatial vector code segments 
transmitted in each of q = I frequency bins where | c <3 > | 

is the length of the code sequence. In this embodiment, an 
information signal b(k) is converted into a code sequence c<:>> 
by a one dimensional trellis encoder, as shown in FIG. 8. Code 



19 



WO 98/09395 



PCT/US97/15384 



segments of length Mt form a spatial vector code c n for a 
single MOFDM frequency bin indexed by n. After training 
symbols are injected, frequency domain interleaving is 
performed by an inter leaver in order to distribute consecutive 
5 spatial vector code segments among well separated frequency 
bins. Interleaving allows the system to exploit the frequency 
diversity of the channel while the spatial coding is a form of 
spatial diversity. 

10 Each of the M T symbols in a given spatial vector code segment 
for a given frequency bin are transmitted from one of the 
antennas. At the receiver, a matched frequency de- inter leaver 
transforms the space- frequency sequence back into a serial 
signal stream. A tilde, ~, above a variable is used to denote 

15 the signal sequence before interleaving and after de- 
interleaving operations. Define 

e (j) _ r <j> -it 
c^' - l c 0 , . . . , c qMT _ 1 J 

20 as the trellis encoder symbol sequence codeword of length qM T 
indexed by j . Further define 

% (q) (k) = [ % ?MR (k) *?M R+ gM T -l(k) ]T 

25 as the received de-interleaved signal sequence due to the 
transmitted code c^) where fM R is the beginning index for the 
received sequence of length qM R spanning q space- frequency 
subchannels. The output sequence due to codeword c ( 3) can now 
be written as 

30 

ft (q) (k) = ft (q) (k)c(3) + ft (q) <k) 

where 
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ft (q) (k) = 



0 ft f+q {k)y 



and 




P T 



{ Z?M T (k) ^ 

\2f +qM T -l < k > / 



The additive noise term % (q) (k) is still white after the 
MOFDM channel operations , 

For a given spatio-temporal symbol code set 
C = { c<l>, c<J)} 

the maximum likelihood detector is given by 
c = arg max P (c U ) I X <q) (k) ) . 

e( D) 



FIG. 7B shows such a detector which is used to generate the 
recovered information stream, £>(k) . Given that the receiver 
noise present in each space- frequency sub-channel is 
multivariate AWGN, it is known that the equivalent decoder 
optimization is 



arg mm 
c(jj 



\-\Jf* «""»(k)o<JI-*<«>(k) = 



This equation can be solved efficiently using a vector Viterbi 
detector similar to that used in ISI channels. The main 
difference here is that while there is correlation in the 
received spatial code segment in each frequency bin, the 
information across frequency bins is uncorrelated . This 
allows the metric computation to be pruned back to the number 
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of states in the trellis encoder at the beginning of each new 
spatial code segment hypothesis test. It is undesirable for 
the encoder to posses parallel transitions because this 
reduces the diversity order of the code to one. Therefore, all 
5 of the encoder input bits are fed to the convolutional encoder 
with rate r and there are there is only one member in each of 
the cosets. 

The inventors have discovered that, given a random Rayleigh 

10 channel process with uncorrelated spatial fading and perfect 
frequency domain interleaving and any code set C, the upper 
bound on the average bit error rate for a 1 x m SIMO channel 
is larger than the bound for a M x mM MIMO channel, even 
though the later transmits data at M times the rate of the 

15 former. This remarkable fact reveals some very interesting 

behavior for the proposed MIMO channel coding structure. 
Although the data rate for the M R = M = M T MIMO channel goes 
up linearly with M, the probability of error bound is smaller 
than that for the SISO channel. While the transmitter power 

20 for each spatial symbol must be reduced as the number of 
antennas and spatial sequences are increased, the length of 
the spatial code segment error vector a n (jl, j2) also 
increases to offset the transmitter power reduction. In 
addition, while the frequency diversity due to the number of 

25 frequency bins spanned by the code error sequence is reduced 
as M increases, the denominator exponent increases due to 
spatial diversity. Thus, as M increases, the effects of 
frequency diversity are replaced by spatial diversity. 
Furthermore, it is clear that the MIMO system can benefit from 

30 additional spatial diversity by setting M R > M T . The m-order 
spatial diversity error performance of a 1 x m SIMO channel 
can be achieved with an M x mM MIMO channel which will again 
achieve M times the data rate of the SIMO channel while 
maintaining lower error probability. 

35 

An alternative code design metric for the spatio-temporal 
coding structure presented in relation to FIGS. 7A and 7B is 
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suggested by observing that the correct error sequence metric 
for code design is clearly the product of Euclidean distances 
for each of the M T length spatial error vector segments in the 
error event sequences. This metric is strikingly similar to 
periodic product distance metrics that are known from other 
contexts . 

An important aspect of the present invention is channel 
estimation. In fast fading channels, overhead penalties for 
conventional multi carrier training techniques can be severe. 
A large number of sub-channels N is desired so that cyclic 
prefix overhead is minimized. Large N corresponds to long 
OFDM symbol duration. Long symbol duration, in turn, requires 
short intervals between training. In conventional channel 
training procedures, an entire OFDM symbol is dedicated co 
training, and several data symbols are inserted between 
training symbols. Thus, a trade-off exists between cyclic 
prefix and training overhead with conventional channel 
estimation techniques . 

Furthermore, in burst-mode transmission applications such as 
wireless ATM, if the average data rate for a virtual circuit 
is low, then the time between ATM packets can be large. In 
such cases, it is not feasible to use an entire DMT symbol for 
training since the channel can change substantially between 
training symbols. What is needed is a training strategy that 
allows "instantaneous updates" for the channel estimation 
algorithm. The present inventors have developed a training 
approach and channel estimation algorithm which injects 
training information along with data into each OFDM symbol. 
The channel estimation algorithm exploits the correlation 
properties of the time varying wireless channel to estimate 
the spatial channel for each MOFDM frequency domain sub- 
channel . 

Imperfect channel knowledge can have an impact on error 
probability, and channel estimation noise in the receiver will 
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limit the performance of a spatio-temporal coding system. The 
inventors have discovered that the effect of channel 
estimation errors can be modeled as an increase in the 
effective noise variance. This noise variance increase is an 
5 interesting function of the time varying channel correlation 
function, the portable velocity, the average channel SNR and 
the design of the channel estimation algorithm. Thus, proper 
design of the channel estimator is critical for low error 
probability communication . 

10 

In all that follows, we again invoke the spatially 
uncorrelated Rayleigh fading condition. Although the spatial 
fading is uncorrelated, there is correlation in the OFDM 
frequency and time domains which we wish to exploit. The 

15 correlation in the frequency domain arises from the delay 
limited nature of the channel impulse response. The 
correlation in the time domain fading arises from the band 
limited Doppler shifts experienced by physical objects which 
move in the vicinity of the portable. We desire a channel 

20 estimation algorithm that exploits these correlation 
properties in an optimal manner. 

To estimate the matrix channel that exists at a given OFDM 
frequency index, note that we can simply estimate the M T 

25 column vectors of dimension 1 x M R . Given that the column 
vectors are assumed to fade independently, an optimal training 
strategy is to transmit M T different training sequences from 
each transmitter antenna and estimate the resulting column 
vectors without considering the information received during 

30 training from the other transmitter antennas. In addition, the 
uncorrelated spatial fading assumption allows each of the 
scalar elements in a given channel column vector to be 
estimated independently. Thus, with M T training sequences 
transmitted independently from each antenna, we can estimate 

35 M R independent frequency domain scalar channel entries. Thus, 
the focus is on a SISO training strategy for the frequency 
domain sub-channels that exist between one transmitter antenna 
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and one receiver antenna. The SISO estimation algorithm is 
directly generalized to the MIMO case by stacking SISO 
estimates from each receiver antenna into columns and 
exploiting the cyclic shift properties of the DFT. 

Our SISO channel estimation strategy will be to transmit 
training symbols in several equally spaced OFDM sub-channels 
with data embedded between training symbols. For a discrete 
time channel which is delay limited to v+1 finite impulse 
response terms, v+1 OFDM training sub-channels are sufficient 
to construct an estimate of all N sub-channels. 

A SISO OFDM channel estimation algorithm is now described. The 
channel frequency domain training symbol sequence is defined 
as 

Z T = diag Z 0 r Z N , Z VN . 

L v+l v+l J 

By construction, 2 T Z^ = Pt^v+1 • 



The channel estimation procedure is as follows. 

1. Given ni past measurements, the present measurement, and n 2 
future measurements of the frequency-domain training sub- 
channel outputs, form ni+n 2 +l measurements of the time varying 
channel impulse response vector ii(v + i) (k) by dividing the 
received known training symbols into the outputs and then 
performing the I DFT operation, i.e. 

£(v+i) (k) = V^ViZt 1 ^ 1 (k) , 
where X$ +1 { k } = \% 0 ( k ) , ( k) *_vN_(k! 

L v+l v+l 



and Fv+i is the 



v+1 point DFT matrix. 
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2. Form the channel impulse response estimate h(v+l) (k) by 
applying an optimal linear MMSE estimation filter 
independently to each of the impulse response measurements 
h< v+1 > (k) , i.e. 

5 

h v+i (k - ni r 

h v+1 (k - n 2 ) 

where ® denotes the Kronecker product and Wh is the scalar 
Weiner filter for h< v+1 > (k) . 

10 

3. Form the complete OFDM channel estimate by zero-padding the 
channel impulse response estimate and performing an N-point 
FFT, i.e. 

15 ft (N) = F N [ h( v+1 ) (k) , 0 0 ]T. 

Given the iid fading assumption on the channel impulse 
response terms, the above channel estimation algorithm is 
optimal within the class of linear MMSE estimators. 

20 

To extend the preceding scalar channel analysis methods to 
estimate the OFDM matrix subchannels, the following procedure 
is employed. Rather than transmitting v+1 training symbols 
spaced by sub-channels, we transmit V+1 frequency domain 

25 sequences, each of length M T . This training scheme is 
illustrated in Table 1 where the notation T(n) represents 
training symbol n and D(k) represents data symbol k. 



h v+l (k) = [ w ; ® i v+1 ] 
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FFT bin 


Content 


0 


T(0) 


: 




M - 1 


T (M-l) 


M 


D(0) 


: 


. 


N/ (V+l) - 1 


D(N/ (V+l) -M-l) 




: 


VN/ (V+l) 


T(VM) 






VN/ (V+l) + M - 1 


T(M(V+1) - 1) 


VN/ (V + l) + M 


D( (V-l) (N/ (V + l) -M) ) 






N - 1 


D(VN/ (V+l) -VM - 1) 



In each of the M T long training sequences, the first symbol is 
transmitted from the first antenna, the second symbol from the 
second antenna, and so on. The first column in the frequency 
domain sub-channel matrix response is then estimated by 
performing the scalar channel estimation algorithm on each of 
the M R antenna outputs associated with the v+l subchannels 
which appear first in the M T long training sequences, and 
stacking the scalar estimates into a vector. The other columns 
of the matrix channel are estimated in a similar manner, with 
the exception that the final frequency domain estimates ft (k) 
obtained from the channel estimation algorithm are cyclic- 
shifted to account for the frequency sub-channel offset for 
each transmitter antenna training sequence. 

Using this technique, for a complex M T x M R DMT channel, only 
M T (v+l) DMT sub-channels are required for training. The 
channel overhead loss due to training and the cyclic prefix is 
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tiien M T (v+l) +v ^ rp^^g training technique only requires FFTs and 

N+v 

FIR channel estimation filters to implement. 

Among the various applications of the present invention, one 
5 of particular utility is a wideband wireless ATM local area 
network for a campus environment. The transmitted digital 
symbol rate is 10 MHz. The portable terminals are mobile with 
a maximum velocity of 70 miles per hour. The RF carrier 
frequency is 5.2 GHz. This application is extremely 

10 challenging for conventional equalizer based communication 
structures due to the large delay spread and extremely high 
Doppler frequency (+/- 540 Hz). The Doppler shift also makes 
conventional CDMA approaches difficult due to the required 
power control loop bandwidth. For these reasons, the 

15 application is an ideal candidate for MOFDM . In one 
embodiment, such a MOFDM system may have 3 transmitter 
antennas and either 3 or 6 receiver antennas. 

Thus, it will be clear to one skilled in the art that the 
2 0 above embodiment may be altered in many ways without departing 
from the scope of the invention. Accordingly, the scope of 
the invention should be determined by the following claims and 
their legal equivalents . 
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CLAIMS 

What is claimed is: 

1. A method of digital wireless communication between a base 
station and a subscriber unit, the method comprising: 
space- frequency encoding a plurality of information signals 
into a sequence of transmitted signal vectors, wherein 
the transmitted signal vectors have M T complex valued 
components and are selected to send the information 
signals over the a collection of independent spatial 
subchannels ; 

transmitting the sequence of transmitted signal vectors over a 

spatial channel coupling an array of M T antenna elements 

at the base station with an array of M R antenna elements 

at the subscriber unit; 
receiving a sequence of received signal vectors at the 

subscriber unit, wherein the received signal vectors have 

M R complex valued components; and 
performing a space- frequency maximum likelihood detection upon 

the received signal vectors to recover the information 

signals . 

2. The method of claim 1 wherein the encoding step comprises 
performing matrix orthogonal frequency division 
multiplexing of the information signals. 

3. The method of claim 1 further comprising the step of 
injecting training sequences into the information 
signals . 

4. The method of claim 1 further comprising adding cyclic 
prefixes to the coded signal prior to the transmitting 
step. 

5. The method of claim 1 wherein the encoding step is 
performed in accordance with a spatio-temporal subchannel 
decomposition of the channel into independent modes. 
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4 

1 6. A digital wireless communication system comprising: 

2 a base station comprising a base station antenna array and a 

3 base station signal processor coupled to the base station 

4 antenna array; 

5 a subscriber unit comprising a subscriber antenna array 

6 coupled through a wireless channel to the base station 

7 antenna array and a subscriber signal processor coupled 

8 to the subscriber antenna array; 

9 wherein the base station signal processor encodes downlink 

10 signal information by matrix orthogonal frequency 

11 division multiplexing the signal information; and 

12 wherein the subscriber signal processor decodes the downlink 

13 signal information by vector maximum likelihood detection 

14 and space- frequency matrix channel estimation. 
15 

1 7. The system of claim 6 wherein the base station signal 

2 processor performs interleaving of the signal information 

3 and wherein the subscriber signal processor performs de- 

4 interleaving. 
5 

1 8. The system of claim 6 wherein the base station signal 

2 processor performs trellis encoding of the signal 

3 information. 
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FIG. 4A 
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